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INTRODDCTION 


The NASA Lewis Research Center has, for the past several 
years, had contracts with Pratt S Whitney Aircraft and 
General Electric to study engines for the Supersonic Cruise 
Airplane Research or SCAR program. Many novel engine 
concepts were considered during these contracts, including 
several that have been broadly termed Variable Cycle Engines 
or VCE's. 

In order to evaluate these new engine concepts and in 
particular as applied to supersonic aircraft, a computer 
code capable of calculating performance of these engines 
throughout the flight envelope was needed. In the past, this 
"matching" of turbofan and turbojet engines was accomplished 
at Lewis with either the GENENG I or GENENG II computer 
codes (refs. 1 and 2). These codes could simulate turbofans 
with up to 3 spools and 3 streams (including aftfans) and 1 
or 2 spool turbojets. It soon became apparent that these two 
codes were not capable of simulating some of the engine 
concepts evolving from the SCAR studies. 

He therefore needed to develop a new computer code in which 
an arbitrary engine configuration consisting of selected 
component combinations could be described at input time. It 
was also necessary to change engine configuration while 
running the code to simulate the operation of various VCE 
concepts, and to optimize the settings of variable 

components such as nozzle or turbine areas (e.g. to minimize 
SFC for a given thrust) . 

Contact with the Naval Air Development Center, Warminster, 
Pa., revealed that they had a computer code, NEPCOMP 
(ref. 3), which already contained some of the features 
desired and whose structure was flexible enough to permit 
the addition of others. This code lacked optimization 
capability and the ability to operate with "stacked" maps 
which would represent variable component performance. 
However, it already had the capability for processing 
arbitrary engine configurations. NASA-Lewis therefore 
contracted with the Naval Air Development Center for the 
joint development of a revised computer code. The objective 
of the joint effort was to develop a code capable of; 
simulating any turbine engine the user could conceive, 
simulating variable component performance, changing airflow 
paths while running, and optimizing variable- geometry 
settings to minimize the specific fuel consumption or 
maximize the thrust. 

An interim version of this new code given the acronym NNEP 
(Navy NASA Engine Program) became operational in May of 1974 
and has been continuously refined since then to include all 
of the desired capabilities. 
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NNEP contains almost all 
the philosophy of constr 
reference 3. The improve 
to NEPCOMP are in th 
optimization capability, 
maps for VCE operation, 
simulate flowpath switch 
configuration schematic 
boattail drag calculat 
format. This present rep 
and provide a summary of 
the code in its present 
its use. 


of the subroutines and incorporates 
uction of NEPCOMP as described in 
ments incorporated in NNEP relative 
e addition of: (1) a performance 

(2) processing of stacked component 

(3) multi-configurations (modes) to 
ing* {^) a computer generated engine 
, (5) throttle dependent inlet and 
ions, and (6) a simpler input data 
ort will discuss these improvements 

the capabilities and limitations of 
form, along with a few examples of 


OPTIMIZATION TECHNIQUE 

As previously mentioned, one of the primary objectives of 
the joint Navy/NASA engine code development was to add the 
capability to optimize the engine performance (e.g. minimize 
SFC for constant thrust). Two basic approaches to the 
optimization problem were investigated: (1) optimization 

inside the loop and (2) optimization outside the loop. 

By "outside the loop" we mean that the engine is first 
matched; then the free variables are changed and the engine 
rematched. This procedure is continued until the 
optimization is achieved. 

By "inside the loop" we mean that at the same time as the 
engine is being matched, the free variables are changing. 
When the match point is finally achieved, the free variables 
will be at their optimum values. Ideally, inside the loop 
optimization should require 2/n times as much computer time 
as outside the loop (where n represents the number of free 
variables). Both methods of optimization were tried with 
results as follows. 

Outside the Loop Optimization 

Five separate methods were tried to evaluate outside the 
loop optimization. These were: 

(1) Hooke-Jeeves pattern search (ref. 4) 

(2) A first-order gradient technique (ref. 5) 

(3) A first-order gradient technique building 

second order information (ref. 5) 

(4) Davidon-Fletcher-Powell penalty function 

method (ref. 6) 

(5) Powell's Principal Axis method (ref. 7) 

The Hooke-Jeeves pattern search failed to find the true 
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optimums. It stopped the search while apparently crossing a 
ridge. The next three methods all require the calculation of 
derivatives by finite difference. Noise in these derivatives 
caused all three methods to fail. The sources of this noise 
are internal convergence loops on thermodynamic properties, 
table lookups, and tolerances on the interface errors within 
which an engine is considered matched. In order to eliminate 
this noise, extremely tight tolerances on convergence loops 
and interface errors would be required and computation time 
would increase significantly. 


Of all the methods tested for outside the loop optimization, 
Powell's Principal Axis method (BOTH) worked the best and is 
the method presently employed in the NNSP computer code. A 
discussion of the computational algorithm used in BOTM is 
given in Appendix A. BOTM is however slow, as probably all 
outside the loop methods will be. Since NNEP itself takes on 
the order of 3 to 7 seconds of CPU time on an IBM 360 
computer to achieve a converged solution and the engine is 
continually rematched while optimizing, computation time 
grows quickly as more free variables are introduced. For the 
two free variable optimization shown in Appendix B, 84 tries 
were required to find the optimum and this required 
approximately 1U0 seconds of CPU time. Since each try is 
near the last converged try NNEP is balancing the engine in 
about 2 seconds per try. As previously mentioned, the 
relatively large amount of computer time required by outside 
of the loop methods prompted the search for an inside the 
loop optimization method. 


Inside the Loop Optimization 

Having successfully incorporated Powell's Principal Axis 
method into NNEP, it was now possible to try inside the loop 
methods and see if they found the optimum which was now 
known. 

Four methods were tried. These were: 

(1) CONMIN (ref. 8) 

(2) Martensson’s method (ref . 9) 

(3) FLEXI - the flexible tolerance method (ref. 10) 

(4) Hamiltonian/ Lagrangian multiplier method (ref. 11) 


CONMIN requires the calculation of gradient information and 
therefore was subject to the same problem of noise as 
outside the loop gradient methods. After consuming much 
computer time without achieving the optimum, the method was 
abandoned. 


Martensson 's 
method with 


method combines the Lagrangian multiplier 
the penalty function method. This method 
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requires the guess of a scalar constant C. Test runs showed 
that for some values of C, equality constraints became 
satisfied but the free variables remained unchanged^ while 
for other values of the free variables 

constraints were not satisfied. It was felt 
would require the determination of its own 
in order to converge. This was deemed 


unacceptable and this method was also abandoned. 


changed but the 
that each engine 
best value of C 
to be totally 


FLEXI does not require the calculation of derivatives. It 
generates a surface of both feasible and near- feasible 
points and proceeds to the optimum by eliminating near- 
feasible points. The near- feasible points are made more 
restrictive untile in the limit, only those points 

satisfying all of the equality and inequality constraints 
are left. In the test problem no progress towards 
convergence was observed. Other researchers in optimization 
theory have noted that FLEXI has great difficulty in 
satisfying equality constraints and therefore no further 
testing of the method was tried. 


The Hamiltonian/ Lagrangian multiplier method was the last 
inside the loop method tried. This method attaches a 
multiplier onto each of the constraint equations essentially 
doubling the number of variables (each control variable will 
have an associated multiplier). The method, however, 
requires the calculation of second partial derivatives which 
are even noisier than the first partials. Optimization 
progressed initially towards the known optimum but as the 
optimum was approached and derivatives became smaller, the 
noise caused the method to fail and no further progress was 
achieved. In addition, calculation of the second partials 
consumed large amounts of computer time. It was therefore 
decided to also drop this method from consideration. 

As a result of the foregoing investigations, Powell*s 
Principal Axis method was adopted for NNEP. 


STACKED MAPS 


Most of the VCE’s evolving from the SCAR studies have to 
some degree variable geometry components ranging from 
variable inlet guide vanes to variable stators and rotors in 
the compressors and turbines. The component maps for these 
variable geometry components represent the component 
performance as a function of the settings of the variable 
features. 


Thus, the 
ratio and 
figure 1 
angle as 


map of corrected airflo 
corrected speed for a 
where there are three 
a fourth parameter. 


w as a function of pressure 
turbine might look like 
separate maps with stator 
NNEP can interrogate this 
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^stacked** map determining the corrected airflow for any 
combination of pressure ratio, corrected speed, and stator 
angle. 

DRAWING OF THE ENGINE 

Subroutine FIGURE was added to the NEPCOWP code to draw a 
schematic of the engine in each of its modes (different- 
airflow paths) when the configuration data is read in. This 
is extremely useful when looking at the code’s output since 
outputs are identified by either flow station number or 
component number. These can thus be referenced to the engine 
schematic previously drawn on the output. 

As can be seen by the example figures shown on the output in 
Appendix B, each time a branch occurs out of the main flow, 
a new column of station numbers and component numbers 
appears. The first component in this new stream is 
identical to the one in the main flowstream where the branch 
took place. The last component is either a nozzle or the 
same component as the one in the main flowstream where 
re-entry took place. 

INSTALLATION EFFECTS 


If desired, inlet and nacelle boattail drag penalties may be 
estimated for the engine, assuming an isolated nacelle, to 
indicate installed engine performance. Inlet drag is 
calculated using combined empirical and theoretical 
relations in which the inlet capture area is sized at the 
design Mach number with a specified inlet bleed requirement. 
At other operating points, the calculated engine demand 
airflow and capture area are used to estimate spillage. 
Inlet spillage drag per unit capture area is then a ssumed to 
be directly proportional to the spillage fraction and a 
full-spillage drag coefficient schedule for the specified 
inlet type. An empirical inlet overboard bleed schedule is 
also used to offset spillage drag by assuming that part of 
the excess captured airflow can recover a fraction of its 
initial momentum. 

Aft-end drag is calculated for the isolated nacelle using 
the mean slope of the boattail section estimated from the 
maximum nacelle diameter and the nozzle exit area setting, 
which varies with power level and airflow throughout the 
flight envelope. An empirical drag coefficient function of 
boattail slope is calculated at each Mach number and can be 
scaled to suit desired aft-end characterist ics. 

Therefore, the installation drag calculations are throttle- 
dependent, require a minimum of inputs, and can be scaled or 
tailored to meet expected characteristics for specific inlet 
types and boattail shapes. 
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PROGRAM DESCRIPTION 


NNEP contains almost all of the subroutines 
the philosophy of construction of NEPCOrtP 
reference 3. This philosophy resulted in 
broken into finite blocks so that the 
desired^ replace individual subroutines with 


and incorporates 
as described in 
a code that was 
user couldy if 
ones of his own 


choosing* The flow diagram for NNEP is shown in figure 2. 

Components 

The individual component types are represented as individual 
subroutines. Engine components fall into two broad 
catagories in addition to controls used to balance the 
engine and optimization variables. 

Flow components - falling under this classification are 


(1) 

in lets 

(2) 

ducts/burners 

(3) 

compressors 

(4) 

turbines 

(5) 

mixers 

(6) 

heat exchangers 

(7) 

splitters 

(8) 

nozzles 

Mechanica 1 components- 

subroutines 

(1) 

shafts 

(2) 

loads 

There is 

a limit of a ti 

of th e 

flow, mechan 


are not represented by separate 


total of 

necha 

variables) allowed within the 
any one type of flow or mechai 
maximum number of controls + o[ 

Subroutine Descr ipt ion 

A brief description cf the function of each subroutine is 
given below. Reference to the NNEP flow diagram (fig. 2) 
indicates the interfacing between the various subroutines. 
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VC ENG 


INPRT 


BOTM 


-is the main routine. It decides 
read input, balance the engine, 
to BOTM for optimization. 

-is the optimization subroutine 
user has the option of printing 
of the engine or only the final 


when to write output, 
or turn control over 


and all printing. The 
each try at balancing 
converged case. 


-is the optimization subroutine which uses Powell*s 
Principal Axis method to find the optimum. Once BOTM 


has been called, it takes over as the supervisory 
routine until an optimum has been found at which time 
control is returned to VCENG. 


CALCFX-is used to evaluate the value of 
rainimized or maximized for BOTM. 


the function being 


NEPCAL^ determines the values of the error matrix used to 
balance the engine, determines the new guesses for 
the independent variables, calls INPUT when directed 
to by VCENG, and calls FLGCAL to perform the engine 
cycle calculat ions* 

INPUT ^reads in all of the input data, and writes out the 
configuration information as determined by CONFIG for 
the various modes onto scratch units. It also calls 
the appropriate data back in when modes are switched. 
At the design point, INPUT calls FIGURE. 


FIGURE-when the configuration data is read in at the design 
point for all of the modes, FIGURE schem at i cal ly 
represents the flowpath on the output sheets. 

CON FIG- processes the engine configuration for each mode. The 
flow components are assembled from inlets to nozzles 
as they would appear in the flow stream. The logic 
to be followed in calculating performance is set by 
CONFIG. 


DMINV -is the IBH 363 
routine used 
derivatives used 


double precision matrix inversion 
to invert the matrix of partial 
in the balancing of the engine. 


FLOCAL-sequentia lly 
order to do 
generated by 


calls the components 
cycle calculations based 
CONFIG . 


in 

on 


the correct 
the flowpath 


INLET -performs 
DBURNR-perf orms 
COMPRS-perf orms 
TURBIN-performs 
MIXER -performs 
HE ATXC-per forms 
NOZZLE -performs 
SPLITR- performs 


inlet calcu lations. 

duct, burner, and afterburner calculations. 

compressor calculat ions . 

turbine calc ula tions. 

mixer calculations. 

heat exchanger calculations. 

nozzle calculat ions . 

splitter calculations (bypass engines). 
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THERM -uses built - in cubic spline curve fits for air, 
stoichiometric combustion products, and water vapor 
to calculate gas properties such as; temperature, 
relative pressure, enthalpy, specific heats, and the 
Universal gas constant. 


TREAD -first is called by INPUT to read in all of the maps 
in tabular form which are to be used by any of the 
components. Then, it is called by each of the 
component subroutines to interrogate the tabular data 
previously read in. 


SPLNQI-is a function used to fit cubic splines 
tabular data being interrogated by TREAD, 
to caculate interpolated or extrapolated 
the tables. 


through the 
Tt is used 
values from 


Computer Code Flow 

Returning now to figure 2 we can follow a typical run 
through the NNEP program. 

Design Point 

VCENG calls NEPCAL which in turn calls INPUT. INPUT reads 
all of the maps from TREAD and then reads in the 
configuration and the cycle data for all of the components 
in all of the modes. This data is then processed by CONFIG 
and an engine schematic drawn by FIGURE for each mode. The 
program returns to NEPCAL which then calls FLOCAL to 
calculate engine performance. Control then passes back to 
VCENG which calls INPRT to print out the design case. 


Of f -Design P oi nt 


VCENG calls NEPCAL which calls INPUT. INPUT detects that the 
point being run is not a design point and the program 
returns to NEPCAL. NEPCAL calls FLOCAL to calculate cycle 
performance. FLOCAL checks after the cycle is calculated 
whether or not the engine is "matched”. If not, 
perturbations are made in each of the control variables to 
generate an error matrix. NEPCAL then calls DMINV to invert 
the matrix. NEPCAL then generates new values for the control 
variables and this process is repeated until a balance is 
achieved. Control then passes back to VCENG ant INPRT prints 
the answers. 

Optimization 

The flowpath followed for a case with optimization is 
identical to that of an off-design case. After the engine is 
balanced and control has returned to VCENG, a check is made 
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to see if optimization is desired. If this is the case^ then 
BOTM is called and takes over complete control of the 
program. BOTM acts as a supervisory routine: perturbing the 
optimization variables; calling NSPCAL which rebalances the 
engine; and, predicts new values for the optimization 
variables. When the engine is both balanced and performance 
optimized, control is returned to VCENG which calls INPBT to 
print the answers. 


CONFIGURING AN ENGINE 

Components are connected together through an indexing system 
which requires numeric coding of each component and flow 
station. Each component can have a primary and a secondary 
upstream flow entering and a primary and secondary 
downstream flow leaving. The CONFIG subroutine searches 
through the components from inlets to nozzles and generates 
the correct sequence of component calculations to be 
performed. This information is mass stored on scratch file 
units numbered the same as the mode;i.e. MODE 1 
configuration data is stored on Unit 1. When a particular 
mode is to be run, the information containing the flowpath 
and configuration data is retrieved from the appropriate 
Unit and this information is processed by the FLOCAL 
su brout ine. 

NNEP uses NAMELIST input as opposed to the fixed field input 
used in NEPCOMP. A typical input card specifying the type 
of component and its position in the flow stream is shown in 
figure 3. As an example of a configuration input card, 
consider a compressor (assigned component number 4) with a 
primary upstream flow station 4. The primary downstream 
flow station number is 5 and secondary downstream (bleed 
stream) station number is 6. Then, the KONFIG input card for 
this example would be as follows: 

KONFIG (1 , 4) = «COMP* ^ 4, 0, 5, 6, 

KONFIG is a doubly subscripted array of dimension 5 X 60. 
Each of the 60 possible components has 5 values associated 
with it. The first value is component type, the second and 
third are the primary and secondary upstream flow station 
numbers and the fourth and the fifth are the primary and 
secondary downstream flow station numbers. Since KONFIG is 
doubly subscripted and we are using NAMELIST, we must say 
KONFIG (1,4)= where the 1 lines up the data correctly for the 
component number 4 (the second number). The zero in the 
example KONFIG card indicates that there is no secondary 
upstream flow station for this component. 

Names of the components are coded as •INLT^, *C0MP*, ’DUCT*, 
•TURB’, ’MIXR’, *HTXC’, ’SPLT’, ’NOZZ*, ’LOAD’, *SHFT*, 


9 


*CNTL’^ *OPTV*, On a UNIVAC 1100 series these would be 

4HINLT, 4HCOMP etc. For loads and shafts which are 
mechanical components, there are no flow station numbers. 
The KONFIG card for a load would just have the component 
name but the KONFIG card for a shaft would have all the 
component numbers connected to the shaft instead of flow 
station numbers. The KONFIG cards for controls (•CNTL’) and 
optimization variables (*OPTV’) are discussed later. 


DEFINING CHAR ACTERISTICS OF COMPONENTS 

Each component type has a separate list of inputs required. 
A typical list of the inputs or specifications is shown in 
figure 4. SPEC is a doubly subscripted array of dimension 
15 X 60. Each component (of which there may be up to 60) has 
up to 1 5 required inputs. Representation of a compressor map 
requires 3 input elements: pressure ratio versus 
corrected airflow versus ”R", and efficiency versus **R»’ 
where ”R*’ represents lines drawn on a typical compressor map 
which roughly parallel the surge line. The introduction of 
the intermediate variable **R** in the map representation was 
necessary to circumvent difficulties in reading the maps in 
regions where two values of corrected flow are possible at a 
given value of pressure ratio and speed. On each map are 
constant corrected speed lines and in addition there may be 
a third dimensional variable if the compressor has variable 
geometry such as stator angle. Each map is given an 
arbitrary unique table reference number so that the computer 
code will know where to look up the map data. 

For a compressor at its design point, the elements of the 
spec array are as follows: (1) represents the value of the 

”R” line passing through the design point, (2) is the 
fraction of the total flow entering the compressor which 
leaves by way of the secondary downstream flowpath (bleed 
flow), (3) (5) (7) and (9) are scale factors which are 
internally calculated by the code to make the values at the 
design point on the map equal the design values for the 
engine being simulated. They should initially be set to 1.0, 
(4) is the map reference number of corrected airflow as a 
function of **R**, speed, and stator angle, (6) is the map of 
efficiency versus "R**, and (8) pressure ratio versus **R”. 
(10) is the value of the stator angle setting, (11) 

represents the fractional horsepower lost when part of the 
bleed is extracted from the middle stages of the compressor, 
(12) and (13) are the desired values of efficiency and 

pressure ratio at the design point on the maps and (14) 
represents the design point value of corrected speed at the 

actual design point on the maps. (15) is not used for 

compressors. 


10 


CONTROLS 


Once an engine has been configured and the necessary 
component information supplied, design point calculations 
may be made to establish appropriate map scale factors. At 
all conditions throughout the operating envelope of the 
engine, flow continuity and an energy balance must exist 
amongst the various components. Those components connected 
by shafts and gearboxes roust rotate in a distinct speed 
relationship. In order to ’’match** the engine at any other 
than design conditions, it is therefore necessary to input 
to the code those component variables that are free to 
change in order to achieve equilibrium. 


This is accomplished through the use of components known as 
’’CONTROLS”. As previously mentioned, a total of 20 CONTROL 
and OPTIMIZATION components are allowed in an engine. A 
typical CONTROL is shown in figure 5. In figure 5 a KONFIG 
card identifies component 30 as a control. There are no 
station numbers for controls. A new input SPCNTL of 
dimension 9 X 60 describes this control. This card is read 
as follows: 


Vary SPEC (1st input) of component (2nd input) so that 
Station Property sub (4th input) at flow station (5th input) 
has a value of (6th input) within a tolerance of +/- (7th 
input). The minimum allowable value of SPEC (1st input) is 
(8th input) and the maximum allowable value is (9th input). 

The 3rd input can be *STAP* for a flow station property, 
’DOUT’ for an output of a component such as static pressure 
difference in a mixer, and ’PERF’ for a performance property 
such as thrust. The meaning of the 4th and 5th inputs 
change as a function of the 3rd input. 

For the case shown here, we will vary the ”R” value on the 
maps for compressor 4 to drive the relative difference 
between the corrected airflow at flow station 10 and the 
amount of corrected airflow that the component downstream of 
station 10 will pass, to zero. Since pressure ratio, 
corrected airflow, and efficiency are all functions of ”R” 
for a compressor, changing ”R” will change all three 
quantities and this will be used to balance the engine. 


OPTIMIZATION VARIABLES 

The last type of component is an optimization variable. As 
shown in figure 6, a KONFIG card for these variables uses 
only the first and fourth positions. Input (1) identifies 
this component as an optimization variable (’OPTV*). Inputs 
(2) and (3) are zero and input (4) indicates which component 
has the free variable. In this example, component 12 has the 


free variable. The SPEC card uses inputs (2) and (3) to 
state the minimum and ma^ciraum allowable values of the free 
variable; input (4) identifies which is free to vary which 
in the example is SPEC(1) of component (12). 


SIMULATION OP TYPICAL VCE 

In this section of the text, a typical application of NNEP 
is illustrated with the simulation of a VCE. 

Figure 7 shows the configuration schematic used to represent 
the engine. In MODE 1, the flows are mixed and exited 
through a single nozzle. In MODE 2, the main and bypass 
flows have been separated, the mixer has been eliminated, 
and a second nozzle has been added downstream of the bypass 
duct. 

As can be seen from this figure, components 1, 2, 3, 4, 5, 
6, 7, 9, 11, 12, 13, and 14 are common to the two modes. 
Component 8, the mixer, is present only in MODE 1. Component 
10 in MODE 1 and component 25 in MODE 2 are the same nozzle. 
The area of this nozzle must change when modes are switched 
as a result of the difference in airflow. This is 
accomplished by using different component numbers in each 
mode indicating a "different” nozzle. Hence, the appropriate 
nozzle area will automatically be used when modes are 
switched. Nozzle 24 is an additional nozzle necessary for 
MODE 2. 

A typical computer output for this engine is shown in 
Appendix B. The carriage control has been turned off to 
compress the output. Input card images appear on the output. 
The first set of input tells how many modes there are and 
which one is the design mode. Next appears a table telling 
which maps have been loaded and how much storage they 
occupy. This is followed by the KONFIG and SPEC cards for 
MODE 1, the computer drawn engine schematic for MODE 1, a 
table of the configuration data and control information for 
MODE 1 and a summary of the input data. 

The KONFIG and SPEC cards for MODE 2, the computer drawing 
of the schematic for MODE 2, the table of the configuration 
data and control information for MODE 2 and a summary of the 
MODE 2 input data follow. 

Since there are only two modes, the code is now ready to 
calculate the design point performance in the design mode 
(MODE 1). A table of updated INPUTS is printed and then the 
design point output. The inputs required to calculate 
installation effects on ail the cases have been turned off 
and therefore installed and uninstalled performance will 
always be the same. 
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The next case calculates the performance for MODE 2. Since 
the only change has been the separation of the flows^ the 
engine is already completely balanced except for the nozzle 
flow. But, we have two new nozzles which will now be 
designed to pass the flow coming into them. The engine thus 
balances without having to iterate. The third case shows 
the performance at Mach 0.8, 36089 feet (11000 m.) at a 
turbine inlet temperature of 2600 (1440 K). For the 
fourth case, we turn on control 29 which varies the TIT so 
that the thrust is 1400 lbs. (5800 n). For the fifth (last) 
case, we leave on the control on thrust and vary the two 
nozzle areas to minimize the specific fuel consumption. Area 
of a nozzle is DAT0UT5 and these have been circled on cases 
four and five to show how the areas were opened to lower the 
SFC (the main nozzle increased by 28 percent and the duct 
nozzle 7 percent) and the SFC has been reduced by over 8 
percen t. 


CONCLUSIONS 

The Navy-NASA Engine Program NNEP has proven itself to be a 
powerful computer code. It can be used to simulate any 
turbine engine made up of combinations of inlets, ducts/ 
burners, compressors, turbines, mixers, heat exchangers, 
splitters, nozzles, shafts, and loads. It can switch modes 
and uses stacked maps to simulate variable cycle engines 
with variable geometry. It has the ability to optimize 
engine performance. The optimization method presently being 
used, however, has been found to be slow and any future code 
development work will be directed at speeding up the 
optimization process or developing a method of predicting 
optimum performance. Additional work is also anticipated in 
installation effects modeling. 

At the present time, the distribution of NNEP is RESTRICTED 
TO GOVERNMENT AGENCIES ONLY. 


APPENDIX A 


BOTH 


Outside of the loop optimization is accomplished by a 
subroutine, adapted from ref. 7. BOTM is considered 

suitable for the problem at hand because it does not require 
derivatives to be calculated. It is significantly faster 
than the better known "one-at-a-t ime" method because it 
systematically generates conjugate search directions — eg., 
along the principal axes of ellipsoidal response contours. 
Ref, 7 shows that for the idealized case in which the 
response contours are actually ellipsoids, the true minimum 
will be located in no more than N iterations (where N is the 
dimension of the problem). As each iteration entails N+ 1 
linear searches, the minimum is found after no more than 
N(N+1) linear searches. Since the contours in some 

neighborhood of a minimum are approximately ellipsoidal even 
for a general non-linear problem, BOTH converges very 
rapidly after reaching this neighborhood. In the early 
going, even after a poor initial approximation, BOTH is at 
least as good as alternate non-derivative methods. 

A typical iteration of the computational algorithm is 
described below and illustrated (for ellipsoidal contours) 
in figure 8. 

Let X (5 = an N-dimension al vector defining the best current 
approximation to the minimum. 

Subscript r = dimensional index, 1 ^ r ^ N 

= N linearly independent search directions in the 
N-dimensional space containing Xp. 

Xr = scalar step length along Y^- direction 

Xy. = current value of X following the r-th linear search 
along Yy- direction. 

Initially Xp is chosen arbitrarily and the search directions 
Yy. are taken to be the coordinate directions. A typical 
iteration then proceeds as follows: 

(1) Choose )vr to minimize f (Xy^j + ^yYy) for r = 1,2,...,N 

(2) Replace Y^. by Y^^. for r = 1 , 2 , . . . , N- 1 

(3) Replace Y^ by (X^^-X^) , 

(4) Choose to minimize f (X^+ /v( X^-X^, ) ) and replace X^, by 

X^+X(X^-Xe) . 

Repeat steps (2) through (4) until the minimum is achieved. 


APPENDIX B 


SAMPLE COMPUTER RUN 
Typical Variable Cycle Engine 


L. Fishbach 
E-8606 
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7,8,fl,n,SPFr.( 1,61 = 3.5, 1, 1,3801, 1 ,380 2, 1 , 1 , . , 1 , 
rt,0^9,-0,^PPr( 1,7 >=?, 2,0^1 t 3 803,1 ,3804d^ U O^ IT- 
9, 0 , 1 1 , 0 , SPFC ( 1 , 9 1 =. 03 , 

, 13,0, 14,0,. SPFC d, 11 1=. 03, 

,6,4, 14,0, SPFC ( 1 ,1 2 1 = 5000, 3 <4 ,0,3«1,0, 

, 7, 2^ O^OtSPFC( 1,131=8000,1, 1,0,0 T I ,1,0,0, 

, SPFC( 1 , 141=-200, 

^pCAiT4.44^§4=l,7,^-5;T4^P«~r-8-r-14-rO, 1-, 

, SPCNJTL ( 1 , 16 1 = 1 ,6, 'STAPdfl, B,0,1 , 

, SPCMTL( 1,471=1,4-, • STAPt-,8,7,0,1 ,l-,2.4, 
,SPrMTL(l ,181 = 1,2, 'STAP*, 8,4, 0,1 ,1 ,2.2, 

, SPCMTL ( l-,4Q|=4-rl, *STAP*,a,-?-,0,4^, 

, SPCMTl Cl ,21 1=1 ,12, 'npUT* ,3 ,12,0,1,0,5300, 

=^4 m4T-i-r-8^4-3 , 0,4 ^-0 , 8 500, 

, 14, 0,1 5,0, SPFC (1,24 1=0,, 98 5, 0,0, .985, 1 ,0,0, I, 
,1 4 , 0,4-2,-0, S P FT 4 l-r'2-5 -rO ,0 , «-98-, l-,0-,-0 , 1 , 

,0,0,24, 0, SPFC ( ! ,26 1=0, 0, 100 , I ,4*0 , . I , 


. 9 , 5600,1 1 
94,5200rd 


B3 


-, ?'r ^ JIT-AP-f-^ f ,44^0^4 1 

KnMPfG(lT2R)=»npyv« ,0^0,25.0, 5*.>5PC( l,2B) = 0,Pr Si00,l ,4«=0,, 1, 

K.QhiE.r^-W2g4.ss . ir 4iT4. a-^-Pf-m--U • P.t«E-Lf :4V1 3aCU43^4)+2-Rna^ 

f.F^'O 

1 .__: - 


I 


4vs>r 

Mod^ 2. 




. <TMIT 

i> 





2 







<ZHMR 

-2> 













C %P1 T 


<SPl T 

3> 



4 


n 


<rnHP 

A> 

' r PMP 

4> 

<piirT 

1 J > 



s 


14 


<TU»R 

A> 

-rniJCT 

*5> 

<MP7 7 

?4> 



7 


1 5 




<T IIPR 

6'> 





ft 






cTllpft 

7> 





o 






'"nurr 

9> 





1 1 






<Mfl7 7 

2S> 





12 




URP( 6) 

Awn 

r.nMP( 4) 

AMP 1 n AP( 1 4 ) AMP 


UP'^ ( 7) 

AMP 

rn^-'P! ?i 

AMP 




OSHAFT (n) TS 
(4 THF i^OLinwUJG KFPRFflFA'TS THF f ONF TG>)R AT T mr, MP.nF= 
F T r T [ r T r ^iR F>jr. I r- f p nr mhm rtq a t t hm pipposfs 


« 

cctdPntitAii- 

—MK.IMn 

■* 

NUHRFR 


” 




1 

1 

» 

7 

4 

• 

ft 

7 

•» 

4 

4 


5 

6. 

2 

5.. 


7 

5 

•* 

Q 

? 

« 

1 1 

2 

*» 

1 ? 

1 1 



11 

»* 

14 - 

10— 


IS 

12 

»V 

1 4 

12-- 


17 

1? 

M 

I P 

1 ? 

» 

19 

12 


2 1 

1 2 

•» 

22 

1? 

•» 

24 

Q 


25 

O 


24 

13 

- 

27 

12 

«0 

2P 

1ft 


29 

1? 

O 




ST AT inws 


IIP-STRFAM 


79 CnMpriNFh’T^ 


nnWNSTRFAM 


T vpp 

STAT TOMS 

STAT TPMS 

I ML FT 

1 

0 

2 

0 

CHMPP.FSR 

2 

0 

ft 

Q - 

SPL’ TTFR 

3 

0 

4 

13 

XnHPRFSR 

^ 

.M . 

_ -5... 

.. ._..-6— 

DUCT P 

F 

n 

7 

0 

TURBINF 

7 

6 

a 

0__ 

TIJPBTMF 

8 

0 

9 

0 

DUCT B 

9 

0 

1 1 

. -0 

nufT P 

13 

0 

14 

n 

SHAF.I_ 

6 . — . 

_ . _4 - _ 


-Q_ . 



mi 


CONTROL IMFORMATION 


15 VARY OATINP 

_44> -VARY_aA-TJNR- 

17 VARY HATIMP 

-4-e V AR Y . n A T IM p - 


1 OF COMPONENT 
-I -nF_COMP-ON€m - 
I OF COMPONENT 
4 - nF CO M PONE N T 


7 SO THAT STATP 
-6 Sn THAT -STATP- 
A SO THAT STATP 


19 

-?4- 


2 ? 

-^-7- 


VARY OATINP 

VAR Y-OAT I ND 

VARY OATINP 
V.A RY- 0A44NP- 


?0 

-OCA SE - 


VARY HATINn 
|-OEAiT4 F ICATTO N — 


I hf COMPONENT 

-1 _0F -x:omponfnt - 

1 OF COMPONENT 
_4 OF . COMPONE NT - 
A OF COMPONENT 
,„F-4-C T T C r OALS— F-NG 


I SO that STATP 
12-SO THAT-OATOUT 
13 SO THAT OATOUT 
3- SO THAT STATP 
5 SO THAT PFPPFT 
TNF oEMPMSTPAT 


R OF FLOW station 11 EQUALS 
H OF FLOW .STA^Tin-N --a-i4WALS— 
B OF FLOW STATION 7 EQUALS 


B OF FLOW STATION 2 EQUALS 
B OF rnMomirAJT 12-4^UAL^ O, 
8 OF C.OMPONFNT 13 EQUALS 0, 
B OF FLOW STATION 14 EQUALS 
't EQUALS 0. 130000 0^ 

I ON -.o.jpaasFS- - 


0.00000 

-0^4DOOOO 

0.00000 

0.00000 

00000— 

00000 

Q.,00000 


Cn^*PnNFNT 
A »q ^— TYPE— 

1 INLET O.lOOOOn 03 
2-f-OMPRE-SR 0^1 80000--01 


INPUT PATA 


_0A4:-IMP^- 


SPL ITTFP 
cnMORi;sJt- 


0.500000 00 
0. 110000-01 


5 nuCT B 0.500000-01 

—6 TUILRIAI^ 0^ X ^0000 . 04- 

7 turbine 0.220000 01 

__0_4>UC4— a 0 . 3 00000^04 

11 nUCT B 0.300000-01 

-4.2 SHA FT O . SO Q OOD 0 4 


0ATINP2 

n. 00000 

-4>^4)0000- 
0.00000 
44^40000-04- . 
0.30000D 00 
0^400000 04 
0. 00000 
— 0^00000 - 
0.00000 
- 0 , 100 0 00 04— 


13 SHAFT 

-14. WnAO - 

15 CONTROL 
L4-— C44NI410L- 

17 control 
.1^ 


0.800000 OA 
-0^^0000 -03- 
0.00000 
-0,440000 

0. no 000 

_44^c:oQ..i. 


1 H 

CnNYP.HL 

0, OOOOO 



0* OOOQQ 

22 

CONTROL 

0.00000 

24 

N077L F 

0,00000 

25 

NOZZLE 

0.00000 


nATAAAP 

0, OOOOO 


fMTROL 

, Ja 44 AP 

c:omt->pl 


0. oocon 

-O^OQOOQ- 

O.OLOOO 


0. 1 00000 01 
4>^00000 
0.00000 
0.00000 - 
0. lonoon 01 
...4Qimoa 04 
0.03000 

o.rOooao 

0.00000 
Q.PftSOOO— 00 
0.980000 00 

-0.^440000 

O.OO'^OO 

- o.noooo 

0.00000 


OATINP 3 

n. 00000 

-0.40000 0 0 1 
0.00000 
■ O.lQQQO a - O l— 
0.00000 

o. ioooon ai 
O.lOOOOn 01 
44.00000 
0.00000 

-^4^00044-04- 
o.ioooon 01 
Q.onono 
0.00000 
-0,4)0000 
, '‘40000 01 
. 40014- 04 

. '.ono 

0.530000-04 
0.850000 04 
o.rQonoa — 
0.00000 
OwiOOOOn -03 
0 *00 non 

0.5QQ0 jn 
0.280000 0^ 


D4TINP4 

o.oooon 

0.376100 

o.oonoo 

0.370700 

0.280000 

0.380100 

0.380300 

o.ooooo 

o.onooo 

-©..44)0000 
0.00000 
0.00000 
0.1 noonn 
o.iQooon 

0.1 OOOOO 

-o.iouoan 
0 . I oooon 
0.400000 

o.ioooon 
- 0.00000 
0.00000 
o.-LQoona. 
Or lOOOOO 
o.l *00-00 
0 . A oooon 


OAT4NP5 

o.onooo 

0.100000 04 
0.00000 
- 0 . 1 nooo-n-04 - 
0.9P0O0n 00 
O. lono0f4--04 

0. loooon 01 

— 0.00000 

o.oooon 
04 — 44.00000 - 
0.00000 

0,00000 

01 0.00000 
o.noooo 

0. OOOOO 
0.00000 . - 
o.oonoo 

04 0. OOOOO 

01 0.00000 

0.r9B5000 00 

0. 980000 00 

04 0.00000 

01 0.00000 
01 0-00000 
01 0. 130000 04^ 


-04 

04- 

OA 

04 

OA 


01 

01 

oi 

01 


OAT 4NP6 
0.980000 00 
0.374200 04 
0.00000 

0.183000 05 
0.380200-04 
0.380400 04 
0,00000 
0.00000 

0.130000 01 

Q.OOOOO 

0.800000 01 
0,000000 441 
0.800000 01 
(4.003000 -04 
0.800000 01 
0.00000 
0.00000 
-0.400000-04 
0.100000 ni 

- 0.00000 

o.eoooon oi 
o.ooooo - - - 
0.00000 


-OATINP7- 
O.OOOOO 
0.100000 01 
0.00000 

o.ooooo 
0.4 00000 01 
O.IOOOOD 01 

0.00000 

0.00000 

o.ioooon 01 
0.00000 
0.00000 
0.00000 
0.00000 
0.00 00 0 - — 
o.onooo 


_ DA T I N P 8 - 

— OAT4NP9 


0.00000 

0. OOOOO 

»• 

0* 376300-04 

0.400000-01^ 


0.00000 

0. OOOOO 


0. 37090 Q.-04 - 

0_ lOOQOO 0 1 

r« 

0.00000 

0. OOOOO 


0.100000 -01 

0.000000-00— 

" 

o.ioooon 01 

0.00000 

'* 

- - 0.00000 — 

0* OOOOO 


0.00000 

0.00000 


Q \ OOOOO 01 

0 oowo 

*• 

0.00000 

0. OOOOO 


0.-00000 

0.00000 


0.00000 

O.IOOOOD 01 


0.00000 

0, iOOOOO OA 



0.00000 

-0.00000 

0.00000 


0.800000 01 0.00000 

0.00000 - 

o.noooo 

— 0.00000 

0.00000 
0.00000 
o.ioooon 01 


0.00000 
0 , 0 0000 
0. OOOOO 
-0.4)0000 
0.00000 


O.IOOOOD 01 

o.ioooon 01 

-Ol- 

0. lOooon 01 

o.ioooon 01 
Qj, i OOP on 00 
o.ioooon 01 
-0.-400000-00- 
0.00000 


THF MAXIMUM COMPONENT NUMBER USED 29 
OW OPF 1 N O W OFI N C USED 


Now 


OOFS nht EQUAL 


25 THF MUMPFP OF COMPONENTS CONFIGUPFO IN ANY ONE MODE - WARNING ONLY 


run/ 


oles 


UPOATEO INPUT DATA TO REFLECT CALCULATEO IN«>UT 
_COMOQNFA4T 


13*^ mpcic - Hode -Plow »y u>cle.^ 


NO. 

— 1 — 


TYPE 

INLFT - 


2 COMPRFSR 
-a-Sm-TTTT^R- 


0ATINP3 
-0*446^60-0?- 
0.180000 01 0.00000 0.000030 06 

^-500000 -00 — O. OOOOO 0*00000 


OATINPl 

0,1000 04T-04 


OATINP? 
- 0 .-0 00 00 


OAT I NP6 

-44.0-0000 

0.376100 06 

0*00000 - 


0ATINP5 

._0,^(V)AO 

0.102390 03 
0*00000— 


OAT TMP6 

0.376200 06 
0.00000 


OATINP? 


0.106620 01 
0,00000 


0ATTNP8 

- 0.00000 

0.376300 06 

-0^00000- 


nATINP9 
-44.4>04400 


0.533650 00 
0,0000(5 


6 COMPPESR 0 . noonn Ol 0 . 36000 n -01 0.66869n 06 0,370700 06 0.372190 02 0.370800 04 0.897910 00 0.370900 06 0.395970 00 


riiJC-T-e a^-^joooo IV- n i n .-a^Aio-ar 

6 TUPRIWF G.35000P 01 0, lOOOOn 

^ — i4jaa-u^i^ — o^2^aoap-4^i — o»^jooo~ 


_ 0 O-- 0^^000 

01 0«38«Jt->n0 00 


o^-ooooo--o4 — o.^rmoo -oo — i 

0.380100 04 0.12005D 01 0.380200 04 


-02 a^ooroo 

0.10141D 01 -0.292V50 00 


8 MIKFR 0.3191^0 

— Q — mjCT - B 0^ 300000- 

10 M0Z7LF n.22l4rP 

-11 DUCT— a Q.aQOQon= 

12 SHAFT o.soooon 

-1-3 SHAEJ. 


14 LOAD -0.200000 
T-IDM - 


03 0.336390 

TU 0^000 OT) - 

03 0.98000P 
.QI_._XUOOOOO- 

04 O.IOOOOO 

03 0.00000 

-F_ICTlCinUS 


04 0.300000 00 0.800000 00 

-0.00000 0..00000 

00 0.00000 0.00000 

0.00000 -— 0.00000 

01 o.ioooon 01 o.ioooon oi 

_a.i_ooooo-Oi 0^00000 0.00000— 

0.00000 0.00000 0.00000 

EMGIME FOR OFMOMSTRATinW PURPOSFS - 


0.00000 0.00000 

o.4Joooa 0..410000 

o.')8ooon 00 o.ioooon 01 

-0^00000 0^00000 

0.00000 o.ioooon 01 


0.00000 

-a.Tioooo- 

0.00000 

_a^Doaa 

o.ioooon 01 


0.00000 

— o.x>oooo 

0.00000 
- 0 . 0000 a 

O.IOOOOO 01 




O.SOOOCM3 00 
-0.410000-. 

O. 000(10 

a.T>aooa- 


-.-Q^ i -oo-onD-ai — OH.-i.4)ooaD~o.i— -a.,aoooQ 

0.00000 0.00000 o.noooo 


O.IOOOOO 01 

- 0 . 00-000 

0.00000 

o^OnfiQO 

0 . oonoo 





V-'omT ( K 

it 





FLOW 

WEIGHT 


TOTAL 

" 

PI nu 


prF?vSURF 



STATPl 


STATP? 


1 

0.-1 000 on 

03 

0.146960 0? 


2 

0 . loooon 

03 

0.144020 02 


3 

onoon _03 

0.288040 02 

r: 

4 

0.666670 

02 

0.2 88 040 02 



0.64 26 70 

02 -- 

— O.l 18100 03 

* 

6 

0. 240000 

01 

0 . 1181 on 03 


7 

0.66179 0 

02 

-0.9 546 50 0? - 


8 

0.685790 

02 

0.551230 02 

* 

.9 — 

88^9-0- 

02 

-- (U 409300- -02 


>iA 


STATTOK' PpnPFRTV OUTPUT DATA 


10 0. 1 0101 o 03 

„ l I 0 . 1 0191 n 03 

12 O.lOlOin 03 

_13 a.- 3-333-3 0 -07 

1 4 0. 333330 0? 


0.394640 02 
0.387800 02 
0 . 3 B? 80 n 02 

0.288040 02 
0.279400 02 


totm. 

TPMpPRATURF 
STATP3 
0.518f.7D 03 
0.518470 03 

0.644650 03 
0.644650 03 
0.100140 04 
0.100140 04 
0.280000 04 
0.245770 04 
0.231070 -04 
0.181260 04 
0.181760 04 
0.181260 04 
0.644650 03 
0.644650 03 


fufl/atr 

PFFFRRED 

MATH 


STATIC 

TNTFRP^CF CORRFCTFO 

RAT (P 

FL-OW 


MUMBFR 


PRFSSURF 


FLOW FRROR 

STATP4 

STATP5 


STATP6 


STATP7 


STATP8 

0.00000 

0.999980 

0? 

0.00000 


n. 00000 


-Q^oonon 

0 -ononn 

0. 102040 

03 

0. 00000 


0. 00000 


- 0.00000 

0.00000 - . 

. -Q.^56A7-90- 

02- 

—0.^0000- 


0.00000 


-n^.nnnno . 

0.00000 

0.3791 90 

02 

0.00000 


o.noooo 


-o.oooon 

0.00000 

-0^11 11 20 

02 

a . 000 00 


0.00000 


-- 0^00000 

0 .00000 

0.00000 


0.00000 


0.00000 


0.00000 

0.297540-01 

0.23 6700 

02 

- o.oonno 


o.oonoo 


-n^onnnn 

0 .286830-01 

0.397990 

02 

o.noooo 


0. 00000 


-0.00000 

-0,28683n--01 

0.91 96 70 

02 

o^3oaaon 

00- 

0.38^200 

02^ 

--0.00000^- 

0.1 91220-01 

0.70944.0 

0? 

0.00000 


0. 00000 


-0.00000 

0.191220-01 

0.73 1390 

02 

o.ioooon 

01 

0.206560 

02 

-0.00000 - . . 

0.1 91220-01 

0. 731390 

02 

0 . 1 19690 

01 

0.146960 

02 

o.onooo 

0.00000 

0. 1 89.6 00 

07 

0.00000 


o.oonoo 


-0.00000 

0.00000 

0. lf^‘^4 60 

02 

0 .709490- 

-0? 

0.386700 

02 

0.00000 




— 

— - 






rOMpPMFMT nuTPUT r'ATA 


.X-DMP-P.NFMT- 
MO. TYPF 


0AT0UT2 


DATPUT3 


0AT0UT4 


0ATOUT5 


0ATniJT6 


0ATnilT7 


OATOUTfl 


0ATnUT9 


2 

CnMPRFSR 

-0,428210 

04 

0.800000 04 0.00000 

O.lROOOn 01 

o.aoooon 04 

o.ioooon 01 

0.102390 

03 

0.900000 00 

0. 200000 

01 

> 


<;pi iTTfR 

n. Rooonn 

on 

n^nnnon n^noono 

- 0.0 0000 

0.00000 

. - 0 . 00000 - 

- 0.-00000 - . 


r\, nnr>r >0 

A- AAAAA 


«• 

4 

comprfsr 

-0. 821810 

04 

0.500000 04 0.00000 

0.1 looon 01 

0.44 8490 04 

o.ioooon 01 

0.372190 

02 

0.880000 00 

0.410000 

01 

*• 

- 5- 

-nUCT-B— 

0^149000- 

00 - 

-0.-500000=01 -0.-300000 -OO 

- 0.?9754D=ai- 

- 0.7333-60- 02 

0.^88400 -04- 

0.300000- 

00 - 

0.1 83000-05 — O.ftqOOOD-00 


6 

TURBIMF 

0.841810 

04 

0.500000 04 o.ioooon 01 

0. 350000 01 

0.3842 80 no 

0.560000 04 

0.1 20050 

01 

0.900000 00 

0. 173190 

01 


7 

tupPT me 

0. 428210 

-04 

0,800000 04- O^innnon-Ol 

0.2 20000-01 

- 0.706760-00 

0.520000 04- 

0.6 302-50 

00 

0.910000 -00 - 



8 

M IXFP 

0.319140 

03 

0.336390 04 0.105980 01 

0. 723470 00 

0.676410 03 

0.882340 01 

0.366450 

03 

-0.9 19920-16 

0. 10219D 

01 

« 

0 

niiPT R 

0 ^ onnnn 


n^^nonnn-ni o-nnnnn 

n-,nonnn 

n^oonnn 

n-,nonnn 

n,nonon 


n , 1 R-^npn ar 

A- AAAAA 


m 

10 

N0Z7LF 

0.714810 

04 

0.225670 04 0.260480 01 

0,230270 0^ 

0.221470 03 

0.980000 00 

0.980000 

00 

0.185320 01 

0. 2604 80 

01 

•» 

. 1 1 

niirT R 

0 , onoon 


-0.300000=01 0 .^) 0000 - -- 

o.ooooa 

0 . 00000 - 

- - 0^00000 Q^Qonno 


0.183000 05 - 

A, AAAAA 


«r 

12 

shaft 

o.oooon 


0.500000 04 0.500000 04 

0.500000 04 

0.500000 04 

0 . 00000 

0.00000 


0.00000 

0. 00000 



1 3 

SHA F T 

n . 00000 


-0..80000D 04 —0.800000—04 

0.800000-04 

- 0.-00000 

- 0.00000 

— 0.00000 


0.00000 

A, O00AA 



14 

LOAD 

- 0.200000 

03 

0.500000 04 -0.00000 

0. 00000 

0.00000 

0.00000 

0.00000 


0.00000 

0.00000 


u 

MATM— n nnnn ai TiTiinP; 

= ,. 0 ^ . R pf;nVFRY= 0.-08-00 - 0 1 

TER ATT QMS 

0 DACCPC 







ATRPi nu.fi ft/^pri 


_ nRn«:c TyRii5:-T 

34 48.12 

-FUEL -F4.nw-I LP/-HR ) 


6883.97 




MET 

THRUST 



7148.12 TSFC 


0.9630 

MET THPUST/ATRFLPW 


71.4812 



- 

' TFITAl TMIJ^ 

riHAri 


A-_nA TJIT-AJ RRAI^P CHAPT HP 

0.-OO- 

RflA TTAI 1 nR-AA- 


0,-00 




INSTALLED THRUST 


7148.12 IMSTALLFD 

TSFC 

0.9630 

SPILLAGE 4 

LIP OPAG 


OoOO 



- 

€0 

MnOF=2t 

EFMO 











‘ 


NOW-RP-I MG--USFD- . - - - 

OCASF lOFNTIFlCATinM FICTICIOUS FNGIMF FPr nFMPNSTR AT TON PURP"'SFS 



Case Z - 

Mode. 2 

, 

FLOW 

— .. WF I-GH T 


• 

STATION 

FLOW 

PRFSSURF 

’* 


S4^ATP4 

S4gVT-PJ 


aT 


jstrs 

ST^TIOM -PPOPFRTY -avFTAUT 



total 

TFMPCPATUPF 
SUTAfP^ 


FUFI./ATR 
PATT n 
_s;.T^TT>^ 


RFFfRR-FO 
FLOW 

— STAT4i§- 


_J4ACH 

NUMRFR 

_S41A^A4^- 




1 

0.1 OOOOO 03 

0.166960 

02 

0-518670 

03 

0.00000 

0.999980 

02 

0.00000 



-> 

--0,400000- 03 

0.446020 

0? 

0-518670 

03 

0.00000 

0. 102040 

03 

0,00000 



3 

n. looooo 03 

0.288040 

02 

0.666650 

03 

0.00000 

0. 560790 

02 

0.00000 



-4 

— o,-666xan -02 

0, 288060 

0? 

0.646650 

03 

0.00000 

0.379190 

02 

0,00000 

“ 


5 

0.6^2670 07 

0. 1 18100 

03 

0.100160 

06 

0.90000 

0. 11 ll 29 

02 

0.00000 



6 

— 0.^*44)000 -4)4 - 

-0,4 48100 

0? 

O.l 0n\6O- 

06 

0.00000 

0,0004)0- 


0,4)0000 



7 

0.661790 02 

0. 956650 

0? 

0. 2'*rooo 

06 

0.297560-01 

0.736700 

0? 

0.00000 

'■ 


8- 

-„a.6857-90 02 

0.551230 

0? 

0.265770 

06 

0.286839-01 

n. 39 7990 

A? 

0,00000 

•" 


9 

0.685790 02 

0.600300 

0? 

0. 2 M 02 0 

06 

0 .286830-01 

0. 519670 

02 

0.00000 

. 


44 

-0,685790 07 

0, 397020 

0? 

0.23-020 

06 

0-286830-01 

0, 535749 

0? 

0 . loooon 



1 ? 

0.6P5790 0? 

0. 397020 

0? 

n.2 -102 0 

06 

0.286839-01 

0. 535740 

02 

0. 122569 



. 

14— 

__-0^ 344330 02_ 

-4). 288 040 

02 

0. 64'*65 0 

03 

0.00000 

0,189600- 

A2 

__ou)oaoo 



1 6 

0.333330 0? 

0. 279600 

02 

0 .66/. 65"^ 

03 

0.00000 

0. 1956 59 

02 

0. 1 OOOOO 

I 


15- - 

— 0.433330 02 

0, 2796-00 

0? 

0.66/ 65 0 

03 

0-00000 

0. 195660 

0? 

0.988 000 




• . . ... 





COMPONENT OUTPUT OATA 




COMPONENT 











— AlC 

T-YPE— 

UATOUT-I 

_ _0ATnUT2 

oatoutt 

0 AT0UT6 

DATnuT-5— - 

^ ...0440 

'*■ 

1 

I NI.FT 

0.00000 

0.00000 


0. OOOOO 


0.190000 01 

0 . looooo 

01 

0.00000 

' 


C.O.MPRFSR 

-0,628210 06^- 

0.800000 

06 

0. OOOOO 


0,180009 01 

0,800000 

04 

n.ioooon 


3 

SPL TTTF-R 

0.500000 00 

0. OOOOO 


0.000-43 


0 .90000 

0. OOOOO 


0.00000 


-4— CHM PILE-SP-^CU 82 1 B 1 n 4)4-- 

-0. 500000 

04 

0.009-^3 


0.110009 01 

0.44 8690 

04 

n.ioooon 


5 

DUCT R 

0. 1 A909P 00 

0.500000- 

-01 

0 .30-3: 3 9 

on 

0 .?Q7560-01 

0.723360 

02 

0.688600 


—6- 

T UR-84 -MF- 

— 0,n417aO-(74 

__0,. 54)0000- 

04 

A. l;)3r 00 

01 

0,350000 01 

_0,3842 80 AO 

0 .-560000 


7 

TURP INF 

0.628290 06 

0. BOOOOn 

06 

n. 100009 

01 

0.220000 01 

0. 706769 

00 

0.520000 


— 9u- 

. duct— 8-- 

— 0.00000 

- 0,300000- 

-01 

0 . 00 3 no 


0. OOOOO 

0.00000- 


- 0.00000 


1 1 

nUCT B 

0. OOOOO 

0. 300000- 

-01 

0.000t)0 


0.00000 

0.00000 


0.00000 

“ 

42- 

S-HAF-T- 

-=0.30-31-70-00-- 

- 0. 500000 

04 

- 0.500-300 

04 

0.500009 06 

0. 500000 

04 

- AUlOOOO 


13 

SHAFT 

-0. 127150 00 

0. 800000 

06 

0.B0030O 

06 

0.800009 06 

0.00000 


0.00000 

♦e. 

14 

LOAO 

_n 9(7r»Dnn q 3 

n t;or>non r>A _n nnn{)n 


_ o^ooooo — 

- 0. 00000 


Ou noon/\ . 

•• 

24 

M0Z71 F 

0. 116320 06 

0.11227D 

06 

0.190) 20 

01 

0.113989 06 

0. 578350 

02 

0.985 oon 


-25- 

_MaZ7L-F— 

— a^-5-3-a7O_0<K- 0.^2 59850 

04 

0.270150 

04 

- 0.265150 06 

0.4A3340 A3 

- 0,980000 


MACH= 0.0000 ALTITUDE 

= 0. 

RFrnvFPY= 

0.9POO 0 1 

TFPATI ONS 


1 PASSES 


SIATIC 
PR FSSURF 

S^ATRT- 

0.00000 
— 0,00000 
0.00000 
0.00000 - 
o.ooono 
- 0 , 00000 - 
0.00000 
0 , 00000 
0.00000 
0.? 16050 

0.166960 
- 0 , 00000 - 
n. 1 67500 
0. 1 66961-) 


I MTFREACE-COIUFCTFD 
FLOW F^ROR 

SJPAW* 


- 0.00000 
0.2?667D-16 
- 0.00000 
- 0.261810-15 
- 0.00000 
-0.r00000- 

0. 107370-08 
0.192330-05 
- 0.00000 
-0,00000 
0.00000 

.^0^4^00 

02 - 0.00000 

OP 0.00000 


02 

02 


OAiniW'T - 
0.9B000O 00 
0.102390 03 
0.00000 
0.372190 02 
0.300000 00 
O, 14K)05n 04- 
0.630250 00 
0.00000 
0.00000 
0,00000 
0.00000 
0 . 00000 


0.9B500P 00 

^o,9aoooo -oo 


-OATOUTR 

0.102860 01 
0.900000 00 
0.00000 
0.880000 00 
0.183000 05 
_T), 900000-00- 
0.910000 00 
0.183 000 05 
0.1 83000 05 
-0.3601-50-^)4- 
-0.296940-06 

0 .00000 

0.189310 01 
a.lft3760--Ol- 


O A- T OU T^ 

0.00000 
0.200000-01 
0. 00000 
0.410000 01 
0.990000 00 
■4)^17.3. 1 -S M > 04- 
0.134680 01 

O. 00000— 

0.00000 

0. OOOOO 

0. OOOOO 

-Q.,00000 


0.190120 01 
O,r^7O150- XU- 


ATPFLOW (LB/SFC) 
_AI-FJr_T44RUST 


TOTAL INLET ORAG 
-J.AI s TA 14. F0-4^HRUST 


100.00 GRPS*^ thrust 

6701,89- TSFf 

0.00 TOTM BPAKF SHAFT HP 

-670U89 — INSTAU FO TSFC 


6701.89 
1-.0272 
-0.63 
- 1,4)272 


FUFL flow (LR/HR) 
T44«UST/AlRFLrW 
BOATTATL ORAG 
- SP I L4.-AGF--^ 44 P 4)RAG 


6883.97 

67,4)1^— 

0.00 

0 . 00 - 




en MAC H=-.^ 8 -rA L -T P=36089 , eT A P- . q 6T S P - FG < 4 f 5 ) = 2/^00- ^74417 


OMGOE 


2 NOW 0FING USFO 

|.P4T^T40N FIC^TICIOUS FNG4NF FOR OF NONSTR-ATIOM PUPPnsFS- 


jCase_3 Mode Z aT ^ut>&onic. crjAise 

STATION property OUTPUT DATA 


- X_XT = 2.600 


FLOW 

- STATfON - 


WFIGHT 
— 4 ^= 0 W — 
STATP 1 


TOTAL 

PRFSSUR4^ 

STATP2 


TOTAL 

- TEMPER AT URF 
STATP 3 


FUFL/ATR RFFFPRFn MACH STATIC INTERFACE COWRECTEO 

RAT^IO- - MUMIVFR PRFSSURF FLOW EP410R 

STATP6 STATP5 STATP6 STATP7 STATP8 




I 


* 

1 

_0 3ZST70--O2- 

0* 3792 AO 

01 

0.3901 90 

03- 

0.00000 

0.1 A 5860 

03- 

-0.8000-00 

on 

o.noooo 


— 0.00000 — 




’ 

2 

0.376800 0^ 

0.481 810 

01 

0.AA026O 

03 

0.00000 

0.10588D 

03 

0.00000 

0.00000 


-0.7B833D-04 





3 

■ 0.3 76800—0? 

-.n ujAAna_n2. 

-0,-45767 0- 

-03 - 

Ji.ooaan. . 





n_ QQOOO 


— OOf>9rt 

• 



' 


0.250470 0? 

0. 10160P 

02 

0^557670 

03 

0.00000 

0.3756A0 

02 

0.00000 


0.00000 


0.659760-07 




* 


_a.2AlA5a- 02 

-0^427540 

02- 

0.873370 

03 

0.00000 

0-.1 07700 

02 

- o.oooog- - 


n,onooo 


rr_0 *00000 




’ 

6 

0.901680 00 

0. A275AP 

02 

0.873370 

03 

0.00000 

0.00000 


0.00000 


0.00000 


0.00000 





- 7 

—D.2A8180 02- 

-0.343820 

02 

. 0.260000 

OA . 

0.2 78780-01 

-0.73 75 OD 

02- 

-0.00000 - 


- 0.0000X1— 


-0.232A8D-Q7- 




• 

8 

0.257200 02 

0.199060 

02 

0.227860 

OA 

0.2687AO-01 

0.397970 

02 

0.00000 


0.00000 


0.182380-08 




•0 

-_9 

— Q^-ZS-ZZOD-g? 

-CU14743a-Q2 

-0.2 13880-04- 0.268740-01 

0.570610 07 o-onono 


0,00000 


~0- nannfi 




“ 

11 

0. 757200 0? 

0. 143010 

07 

0.213880 

OA 

0.7687AO-01 

0.536710 

02 

O.IOOOOO 

01 

0.776A10 

01 

-0.7370A0-08 





. 12 

o^zs-izon 02 

0. 1A301D 

02 

0.213880 

OA 

0.2687AD-Q1 

0. 53 67-ID 

02 . 

- Q.LAA65D 

01 

0.3292AD 

01 

O.OOOOQ 




•* 

13 

0.126330 02 

0. 10160D 

0? 

0.557620 

03 

0.00000 

0.1 89A7D 

02 

0.00000 


0.00000 


-0.00000 




’* 

lA 

—0. 12633a 02. 

0. 985520 

01 

0.557620 

03 

0.00000 

0. 19533D 

02 

-O.IOOOOO 

01 

0.5 20A6D. 

0 1 

0.850AB0-Q7 




'* 

IS 

0. 126330 0? 

0. 985S20 

01 

9.557620 

03 

0.00000 

0. 195330 

02 

0.175110 

01 

0. 3292AP 

01 

0.00000 










COMPnNFNT OUTPUT DATA 











COMPONEMT 

















•• 

MO. TYPE 

OATOUTl 

oathut? 

0ATnUT3 

OATOUTA 

06T9UT5 

0AT0UT6 

OATHUT 7 

OATOUT8 

DAT0UT9 



.1 INL£T - 

. 0*90713D 03 

0. 77A60D 

03 

0.A5B95D 

03 

0. 11 2830 01 

0. 152 AAD 

01 

0.300000 

00 

0,960000 

00 

0.97295D 00 

0- 3608 90 

05 



2 COMPRFSR 

-0. 1A9960 OA 

0. 77A77n 

04 

0.00000 


0.181910 01 

0. 800009 

OA 

0. 10532D 

01 

0.1 07309 

03 

0.891A0D 00 

0.2108 70 

01 



3 SPI TTTFR 

0^Sn43fln_ OQ ■ 

JCUQCtOOO-. 


. 0.00000 


0.00000 

... 0.00000 


■—0. 00000 -■ 


-.0.00000 


O.OOOOQ _ 

n. oonnn 




A COMPPESP 

-0.270930 OA 

0.463340 

OA 

0. 00909 


9.10A760 01 

0. AARA90 

OA 

0.996390 

no 

0.372190 

02 

0.888040 00 

0.A208 in 

01 



5-_.nUCT_B- 

-Q*153A9D -00 - 

Q.5QOOOD 

-01 

0.300000 

00 

0.27878D-01 

0.72336D 

02 

0.242320 

OA 

0.2 8 8150 

00 

0.183000 05 

0.99QQQD 

00 



6 TURPI NF 

0.290930 OA 

0. 463340 

OA 

0. 100000 

01 

0.3A8390 01 

O.^PA?80 

00 

0.53B5AO 

OA 

0.1 20050 

01 

0.898090 00 

0. 172720 

01 



7- -TUP.aiM£ 

. -Q.1A995D OA 

0.776270 

CA 

0.100000 

01 

0-22120D 01 

0.706760 

00 

0.52A03D 

OA 

0.63025D 

00 

0.909990 no 

0.I3502D 

QI_ 

>* 


9 DUCT B 

0.00000 

0. 300000 

-01 

0.00000 


0.00000 

0.00000 


0.00000 


0. 00000 


0.183000 05 

0.00000 


> 


11 niiGT P 

n-,nonnn 

- CU30Q00D- 

-ai 

0.00000 


0.00000 

. 0.00000- 


_ O.OOOOQ 


0.00000 


_0. 1 8303 D 05- _ 

-0.00000 




12 SHAFT 

-0. 1 A728O-02 

0.46334^ 

OA 

0.A633A9 

OA 

0. 463340 04 

0. A6 33A9 

OA 

0.00000 


0.00000 


-0.506750-06 

0.00000 


- 

* 

13--- SHAFT-. 

—O, 65 8AS 0^-04- 

0.7762 70 

04 

0.776779 

OA 

0. 776770 04 

o.onooo 


0.00000 


0^00000 


0. A3913D-07 

Q.OOOQQ_ 





lA LPAO 

-0.700000 03 

0.463340 

OA 

0.00000 


n. 00000 

0. 00000 


0.00000 


o.ooono 


o.onooo 

0.00000 



» 

2-A-..-NQZZLE- 

_0. 5J.886D -03- 

0. 132140 

OA 

0.29933D 

01 

0. 1341 50 OA 

0.57R35D 

07 

0.98500D 

00 

0.98500D 

00 

0.189360 01 

0.29 933 D 

QJ 


u 

25 NOZZLE 

0.236270 OA 

0.295500 

OU 

0.434359 

01 

0.301539 OA 

0. 1633 AO 

03 

0.980000 

00 

0.980000 

00 

0. 1 «A190 01 

0. A3A35D 

01 

** 


MACH= n-Aonn ai titiiof 

- 36089. 

R ECPVFRY = 

0. 96 00 13 1 

TERATT DMS 

30 PASSES 









AIR FLOW 11 B/SEGl 

37-. 

68 

GROSS THRUST 

2881. 

-07 

_ FUEL FLOW 

fLR/HRl 


_ 24 73.18 





NET THRUST 


1973. 

89 

TSFC 



1.2276 

MET THRUST/AToflow 


52.3900 



- 

*0 

TOTAl IMl FT 

OQAG 

_ 90.7. 

1 8 

TOTAl RPAJ4E 

.SHA.£I_H8 

-0. 

-OO 

ROATTAII nRAG 


n,no 




- 

TNSTALLFO THRUST 

1973. 

89 

INSTALL EO TSFC 

1.2276 

SPILLAGF 

•- UP DRAG 


0.00 



•< 


“ f.n SPFC(9,29)=:.l,LABFL=T GFNO 

*• OMDQf 2- hlQ^. _aEI -USED 

“ warning EXIT VELOCITY IS SONIC *CnMPONFNT 5 

- WARNING FXTT M FI nO I T Y_ Jf <; <;nNTr, »rn MPnNFNT 5 . 

*’ WARNING EXIT VELOCITY IS SONIC ♦COMPONENT 5 

- WARRING. ♦ t<L-EXiX-V£U3CLTY— IS-SnJ^lC-J^COHPOMEMT 

- warning ♦♦♦ FXTT VELOCITY IS SONIC ♦COMPONENT 5 

“ WARRIAiG--iLt*.-£XIT-V£LDC-l4'Y--TS-SnAlIX.--^CnMPnNFNLT 5 

•’ WARNING ♦♦♦ EXIT VELOCITY IS SONIC ♦COMPONENT F 

“ WARNING AAA FXTT VFLOCITY IS SONIC ACOMRORgRT- — 5 ^ 

WARNING ♦♦♦ EXIT VELOCITY IS SONIC ♦COMPONENT 5 
“ Q C A S-P- 1 D ENT I-F- ir -AX-I-ON VARY— THF T I T TT)-MAKF THF THRUST -1 300 +/- .4 LRS, 

- C ase . M - T’Vji-o'TTTe. _ bac< so "thrwsT = I 3 . 

*’ STATION PROPERTY OUTPUT DATA 


FLOW 

SJ-ATI-ON 


WEIGHT total 

__E4 nw PRESSURE 

STATRl STATP? 


TOTAL EUEL/AIR 

TFMRFTIATURE — RATT43 
STATP3 STATP4 


REFERRED MACH STATIC INTERFACE CORRECTED 

E4.0W number .. - -PRESSURE E-LOW E^kOR 

STATPS STATP6 STATP7 STAT«*F 



Or ; 0- ■ * 


C‘. P . 

. 

■0.2024 


>. 325070 02 

6 

0.754900 

00 

0.325070 0? 

. _7 

-0.^-06530 

02 

0. 26277D- 02 

8 

0.214070 

07 

0.151160 0? 

-9 

0^2.1A07P- 

-02- 

T7OO--02 


11 

-12 

-1-4 


0.?l407n 02 
X), 24 4070 X42 
0.114550 0? 
_0-,4 14-SBn 0? 
0. I 1455P 0? 


0.10B45O 0? 

0,108^50 0? 

0.90325n 01 
0.874150 01 
n.87M50 01 


G • ■ ?■ j ' 'J n 

0«7951%O 03 
0,795130 03 
0.219380- 04 
0.191290 04 
41^174-970 04- 
0.178070 04 
0.178070 04 
0.535360 03 
0.535^40 03 
0. '^35360 03 


O..Cj^‘j0^cr - - - 

0.00000 

-(U- 000 00 

0 .00000 
0,216430-01 
0.20864D-01 
--0.40464-0-01- 
0.208640-01 
-0. 208640-01 
0.00000 
- 0.00000 
0.00000 


.. s. . ; 7i; 0/1 

- XI . 04 — 

0,546040 02 
-0.443220 02 
0.00000 
0.23^540-02- 
0.399700 02 
-O.-^3l40-02- 
0.530360 02 
0.539360 0? 
0.189520 02 
0.1953 80 0 2 
0.195380 0? 


V T fr-GOOOO-OO 

Oe ooono 

- Q ., 40 00 0- 

0.00000 

0.00000 

0.00000 

-0.410000 — 

0.00000 

_a.xioooo 

O.IOOOOD 01 
0,14157D -04 

o.oooon 
0.100000 01 
0.119150 01 


CnMPPMFN'T niiTPLil data 


CO A3 POMP NT - 
NP. TYPP 


OATnUTl 


0ATnuT2 


D/‘TPUT3 


0AT01IT4 


nATPUT5 


0ATnUT6 


- -Cj o ■' .. 

0.00000 

_„f>.000(ic 

0.00000 

- 0.00000 
0.00000 

- 0.0 0000 

0.00000 
— 0.00000 
0.584680 01 
0.329240 01 
0.00000 
0.462680 01 
0.329240 01 


OATPUTT 


-4 c '>0 :^30 
-0.8361 I n-06 

- 0.00000 

-0.324220-06 

-o.oonoo 

0.00000 

-0.182100-04 

-0.541860-06 

-0.00000- 

-0.255150-06 

o.oonoo 

- 0.00000 

0.749060-05 

0.00000 


0ATnilT8 


0AT0UT9 


1- 

. U«LFT -0.7-80710 03 

0.774680 

03 

0.458950 

03 

0.1 12830 

01 

0. 152440 

01 

0.800000 

on 

0.960000 

00 

0.872950 

00 

0.360890 T)5- - 

2 

COMPPESR -0.105630 04 

0.680930 

04 

0.00000 


0.1 3027n 

01 

0 . Roooon 

04 

0.923860 

00 

0.1 02390 

03 

0,901640 

00 

0. 187470 

01 


SPL IT^*^® n PO 

0. QOOOO 


0.00000 


— 0.00000 


-0.00000- 


0.00000 


- 0^00000 . - 


0.00000 


0, OOOOO 


4 

CHmppfsR -0.185990 04 

0.425380 

04 

0.00000 


0.116170 

01 

0. 44 84 90 

04 

0.932 700 

00 

0.3721 90 

02 

0.907140 

00 

0.359890 

01 

5 04JCT-8 0.4A9090-00 

0.^500000- 

-01 

0.300000 

00 

0.216430- 

•01 

0.723060 

02 

0.1575OD 

04 

0. 301740 

00 

0.1 83000 

05 

0, 990000- 

-OO- 

6 

TURBINE 0.205990 04 

0.425380 

04 

0 . loonop 

01 

0.352240 

01 

0.384280 

00 

0.533240 

04 

0. 1 20050 

01 

0.898130 

00 

0. 173840 

01 

' 7- 

__TUR8IN^ -0, 105639-04 - 

-0.680-930 

04 

0, 10000 n 

01 

0.222240 

01 

0.706760 

00 

0.501690 

04 

P.630250 

00 

0.911370 

OO 

0.^1353*20-03 

9 

nuCT P 0.00000 

0.300000- 

-01 

0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


0.183000 

05 

0.00000 


LI 

nnry p ^ ^OOpo 

-Ow-ROOOOO* 

-OL - 

0.-00000 


- 0.00000 


—0.00000 


_ 0^00000 


- 0 , ooooo 


-0.1-83000 

05 

P, 00000 

> 

12 

SHAFT 0.470180-01 

0.425380 

04 

0.425380 

04 

0.47538P 

04 

0.425380 

04 

0.00000 


0 . ooono 


0.228250- 

-04 

0.00000 

• 

-L3- 

— -SHAF7 0-,7XX6390-0L 

0.680930 

04 

0.680930 

04 

- 0.680980 

04 

0.00000 


0.00000 


o.oonoo 


0.19 5390- 

-04 

- 0 . ooooo 

* 

14 

LPAO -0.200000 03 

0.425380 

04 

0.OOP''O 


0.00000 


0_QaQDQ_ 


0.00000 


0.00000 


0.00000 


0.00000 


' -2^ 

_J4077LE 0..4 3947 0-03- 

-O. 123440 

04 

0.2661 1 0 

PI 

0,175320 

04 

57 83 50 

02*^ 

1 0.985000 

00 - 

0.9 8500P 

00 

0.1 89360 
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